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Abstract 

We investigate the effect of long range magnetic interactions on the renormaUzation 
group (RG) evolution of local Cooper pairing interactions near the Fermi surface in high 
density QCD. We use an explicit cut-off on momentum modes, with special emphasis 
on screening effects such as Landau damping, to derive the RG equations in a gauge 
invariant, weak coupling expansion. We obtain the Landau pole A ~ /^^"^ ' 
although the structure of our equations differs from previous results. We also investi- 
gate the gap equation, including condensates of higher angular momentum. We show 
that rotational invariance is unbroken at asymptotically high density, and verify that 
A is the correct value of the gap when higher modes are included in the analysis. 
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1 Introduction 



In this paper we study the behavior of quark matter at high density and low temperature 



[|l|-||T8|. Under these conditions QCD exhibits color superconductivity, caused by the con- 
densation of diquarks. This condensation is analogous to the Cooper pairing observed in 
ordinary superfluids, and can be shown to occur in the presence of even very weakly attrac- 
tive interactions. The reason for this is the special nature of physics very close to the Fermi 
surface (FS). 

At high density, characterized by a chemical potential /i which is much larger than the cur- 
rent quark masses and the QCD scale Aqcd, the typical momentum transfer in quark-quark 
interactions is of order fi, and therefore it is plausible that the dynamics can be understood 
via perturbative gluon exchange. Of course, small angle scatterings, which involve small 
momentum transfers, are still problematic and require special attention. 

In recent work. Son [0 showed that long range magnetic interactions lead to a modifi- 
cation of the RG equations originally derived in [^] for the case of local interactions. The 
magnetic effects are strong enough to modify the parametric dependence of the position of 
the Landau pole (and hence the superfluid gap A) on the gauge coupling constant. Son finds 

A ~/i(7-5exp(--^) . (1) 

That A should scale like exp{—c/g) is easy to see [§] by considering the gap equation with 
a massless, or weakly damped, gauge propagator. In the usual case of a local four fermion 
interaction, the gap integral exhibits a logarithmic divergence which, roughly speaking, is cut 
off near the FS by the gap itself. Solving for A yields a result of the form A ~ exp(— c/G), 
where G is the four fermion coupling, and is of order g"^ if it arises from the exchange of a 
gauge boson. However, if the four fermion interaction is replaced by a weakly damped gauge 
propagator, an additional logarithmic divergence appears due to small angle scattering of 
the fermions. This divergence is again regulated by a scale related to the gap itself, and the 
resulting exponent of the solution is roughly the square root of the what appeared in the 
local case: A ~ exp{—c/g). 

Our intention here is to understand this behavior in terms of the evolution of Cooper 
pairing interactions near the Fermi surface. Despite the long range of the magnetic inter- 
actions, we find that the problem can still be formulated in terms of local operators which 
are, essentially, terms in the expansion of the magnetic gluon form factor. The reason that 
this is possible is because we retain at all times an explicit cutoff on long wavelength modes, 
which keeps all quantities finite. This approach is somewhat different from that of Son |]^, 
who studied the RG evolution of scattering amplitudes themselves. As discussed below. 
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our results differ from his. Primarily, we believe that this is due to conceptual problems in 
applying the RG directly to scattering amplitudes. 

The paper is organized as follows. In section 2 we give a description of our cutoff scheme 
and the resulting effective lagrangian. In section 3 we discuss the problem of small angle 
scattering and gluon screening effects. In section 4 we compute our RG equations and 
compare them with Son's. In section 5 we investigate the possibility of breaking rotational 
symmetry. The final section contains additional discussion of our results. 



2 RG Scheme 



We adopt a Wilsonian RG procedure, with a hard IR cutoff on spatial momentum, A ||5|, p!9 



As A — > 0, only quark excitations very near the FS, as well as soft gluons, are left in the 
effective theory. Our prescription differs from what is often used in QFT, where a hard cut- 
off is imposed on energy as well as momentum, but it has some advantages. In particular, in 
our scheme integration over modes corresponds to shrinking the Hilbert space of the model 
in the basis of energy eigenstates. 

The effective Lagrangian has the form 

U = CqcD + T.^n (2) 

n 

where the On are local operators involving quark and gluon fields, which are the result of the 
integration over higher frequency quarks and gluons. The most important of these operators 
are the marginal Cooper pairing interactions which involve four quarks. All other quark 
interactions can be shown to be irrelevant in the limit of small A 0. 

The Cooper pairing interaction is of the form 

G{k - q)ip+{ko, k)-fi,Pmiip+{qo, q)ip+{ko, -k)YPL,R^+{qo, -q) , (3) 
where ^+ denotes the projection of the quark field 

V'+ = + , (4) 

and consists of quark, rather than antiquark degrees of freedom. We will be interested in the 
case where all of the external quarks in this operator are essentially on-shell. Note that the 
incoming and outgoing quarks have almost equal and opposite momenta. Near the FS, the 
form factor G{k — q) becomes a function of angle 9 = j^, since \k\ ~ |g| ~ fi, and /cq, qo — 0. 
(Strictly speaking, in the case of Landau damping, it is a function of energy as well as angle. 



3 



We will always assume that the energy transfer ko — qq in the gluon line is much less than 
but of order A.) In this paper we restrict ourselves to the 3 color channel, which is attractive 
and to the LL (or equivalently RR) chirality channel, which has been shown to dominate the 
LR channel [^. It is straightforward to derive the related RG equations for other channels 
using our techniques. 

In previous work we took the form factor G{6) to be a constant This is appropriate 
at sufficiently low energy, if screening masses exist for both the magnetic and electric compo- 
nents of the gluon, which is likely to be the case at intermediate densities where the coupling 
is not small. However, as argued by Son 0, a magnetic mass for the gluon is unlikely to arise 
within perturbation theory. The magnetic mass is due solely to nonperturbative effects, and 
is presumably of order exp(— l/gf^). At high density it is therefore likely to be smaller than 
the eventual superfluid gap, and hence plays no role in the analysis. Instead, focused on 
the role of Landau damping on the magnetic interactions. 

In order to consider long range magnetic interactions, it is necessary to expand the form 
factor G{d) in components with definite angular momentum. We can then study the RG 
evolution of each of these components. Let 

G{d) = ^{2l + l)Pi{cosd)Gi , (5) 
I 

so that each component 

Gi = - d{cos9) Pi{cos9)G{9) . (6) 
2 j-i 

This integral exhibits a logarithmic divergence in the case of a massless, or weakly damped, 
gluon (figure 0). However, in our regularization scheme G{6) contains only the effects of 
gluons which have been integrated out above the cutoff A. For nonzero cutoff the compo- 
nents Gi will be finite, but exhibit a logarithmic dependence on A[]. It is this logarithmic 
dependence that leads to the constant Oi^g"^) term in the RG equations noted by Son [0]. 

3 Screening Effects 

In this section we discuss the incorporation of screening effects in our procedure. As men- 
tioned previously, small angle scattering of quarks must be considered carefully, as the simple 
perturbative expansion in powers of g"^ may break down. Indeed, it is easy to see that vac- 
uum polarization corrections to any graph can become important if the gluon momentum is 

^The precise form of this logarithm is dependent on the IR behavior of the propagator, as we wiU see 
below. 
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Figure 1: One gluon exchange. 



sufficiently small. Resummation of these effects leads to a screened propagator which is well 
behaved at small momentum, as long as the energy is nonzero. 

The gluon propagator, including vacuum polarization effects from virtual quarks, 
has the following form in covariant gauge: 



where Q = {q^, q) = (— tu, q) is the gluon Euclidean 4-momentum, and and P^^, are 
transverse and longitudinal projectors correspondingly. In our leading order calculations 
the propagator will always appear contracted with gamma matrices next to on-shell external 
quark hues. Thus the gauge dependent part of (0) will vanish due to the equations of motion, 
leading to a gauge-invariant result. Henceforth we will simply set ^ = in (0). 



The functions G and F are related to the gluon polarization tensor 11 

F — 1144 , 



which reflects the fact that the Lorentz symmetry is broken to 3D rotational symmetry. The 
explicit expression of U^^ to one loop is 

n,.(Q) = 9"J^, Tr[^. ^7. (^- m A{K)A{K - Q) , (8) 



where A{K) = 



Figure 2: Vacuum polarization correction to gluon propagator. 



If the energy and momentum transfer in figure (0) are small then one may neglect Q 
in the numerator of (||), as the dominant momenta in the loop will be ~ /i (this is the 
equivalent of the hard thermal loop approximation). Then F and G take the familiar form 

m 



(1 - t^»(t)) 



g2 V q q 



n 2 
G = rrij^ — 



(l - (-) )^o(-) + - , (9) 



where 



L.{x) = ^ln|^ , (10) 

and m\j = is the Debye screening mass. The small x expansion of G leads to the 

Landau damped magnetic gluon propagator 

while the expansion of F leads to the usual longitudinal propagator, with Debye screening. 
The effect of Landau damping is to cut off the small-q divergence in ([Tl|) at g ~ Qo^^^'d'^- 

One must be careful to compute the screening effects using a Wilsonian cutoff. When 
we integrate over a shell in momentum space only the contribution of quarks which have 
already been integrated out is to be included in the screening effects. This means that we 
must re-examine the calculation which leads to the Debye mass and Landau damping, and 
use cutoff-dependent versions of the vacuum polarization = 11^^, in our RG. 
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The terms in F and G can be shown to result from integration of quark modes within 
roughly q of the FS. For example, the Landau damping term originates from an integral of 
the form (arising from (^ |20]): 



where Ei = k , E2 = \k — q\ 



dk dQ 



E1E2 
k — q cos 9 and 

n{E) 



n{E,) - n{E2) 
iuj — El + E2 



(12) 



1 



f3{E-f,) 



(13) 



is the Fermi-Dirac distribution for quarks. n{E) is a theta function of (/x — E) in the low 
temperature limit in which we work. For a given value of 6 

to Landau damping must be within ~ q cos 9 of the FS due to the theta functions. 



the modes which contribute 



If the momentum transfer q < A << fi, then none of these modes are included in the 
calculation of n^^. The corresponding E^ and are changed drastically: in particular 
the leading /x^ behavior of F vanishes, and G becomes proportional to /x^. Fortunately, our 
interest is only in the logarithmic divergences of diagrams, which are dominated by gluon 
momenta satisfying the limit A << q << fi, due to the form of Landau damping, (i.e., 
A ~ A, while the dominant momentum transfer is g ~ A^/^m^^.) In this limit the results 
for F an G given in (H) are accurate. 

Finally, we mention the issue of screening due to the diquark condensate itself, which is 
necessary for a self-consistent description of the region near the FS [|16|, [l^. In a conventional 
superconductor the relative size of the magnetic penetration depth A and the correlation 
length C, ~ A~^ determine whether one is in the type I (^ >> A) or type II (A << ^) 
regime. In a type II superconductor the magnetic screening length can be computed using 
the London formula, and is proportional to the total density of superconducting particles. In 
our case this would lead to a rather large effective screening mass A~^ ~ g'^fi'^ relative to the 
gap size. However, in high density quark matter we are actually in the type I, or Pippard, 
regime. In this regime the effective screening mass is much smaller, and scales with the gap 
A. A direct calculation of the gluon vacuum polarization using the quark propagator in the 
presence of a gap (see (|3T|) below) shows that the London-type screening applies only to 
long wavelength gluons with momentum less than ~ A. Harder gluons experience a much 
smaller screening of the Pippard type. It is easy to see that these effects are too small to 
affect our RG calculation; their contributions are dominated by Landau damping effects. 
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4 RG Equations 



To obtain the RG equations we need to evaluate the scale-dependent quantum corrections 
to the form-factor G{k — q). Let us consider the effect of integrating out quark states in the 
momentum shell A' < \q\ — fi < A, and gluon states in the momentum shell A' < |g| < A. 
The tree level contribution comes directly from one-gluon exchange, figure |I|, while the 
one loop contribution comes from the box diagram, figure ^ It is important to note that 
only one topology of the box diagram contributes. The diagram with "crossed" gluon lines 
does not have the same kinematic structure as an iterated Cooper pairing interaction, and 
is subleading. Actually, in our effective theory the box diagram contains several different 
contributions. The most important contains two local four fermion interactions, and is 
actually a "bubble" diagram, with running form factor coefficients. These coefficients contain 
the effects of previously integrated shells of quarks and gluons. The other contributions 
involve at least one exchange of a soft gluon within the momentum shell, and are suppressed 
in the thin-shell limit. Thus, as we discuss below, the result of the one loop part of our 
calculation is essentially the same as iterating bubbles with form factor vertices. 

Solving the RG equations is equivalent to summing up an infinite series of ladder graphs 
corresponding to nearly colinear scattering mediated by gluon exchange. This corresponds 
to the "rainbow" approximation in which the gap equation is solved in section 5. 




Figure 3: Box diagram. 

Let us elaborate on why the tree and box graphs are sufficient to compute the leading 
order solution to the RG equations. Consider integrating out all modes between an upper 
cutoff ~ fi and a lower cutoff given by A. For simplicity, we can consider doing this in a 
single step, rather than shell by shell. A systematic expansion for the result is possible in 
powers of the gauge coupling g and powers of t ~ In A. Since the RG evolution terminates 
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at a Landau pole of order ~ exp(— 0/(7), t ~ l/gf in our counting scheme. As we shall see 
below, (figure |^) is of order gH, and so contributes a term of order g'^ to the RG equation. 
The box diagram (figure ^ is of order gH^ and so contributes terms of order gH'^ ~ g'^. 
Higher order loop corrections to either of these graphs (e.g. from vertex or wavefunction 
renormalization) are suppressed by at least gH ~ g. Of course, as discussed in the previous 
section, the screened gluon propagators must be used in our computations, since some of the 
logarithms which arise are due to small angle scattering. Additional radiative corrections 
beyond vacuum polarization effects on the gluon propagators are suppressed by at least a 
power of g^t. 

Because the rotational SO (4) symmetry is broken to SO (3), the coupling in (||) will split 
into its "temporal" G4 and "spatial" Gi parts, which we refer to as A and B components, 
according to the notation previously used in 

A{k - q)^+iko, k)-f4PL^/j+{qo, g)V'+(fco, -k)-f'^PL^+{qo, -q) + 
B{k-q)tlj+{ko,kh,PLtlJ+{qo,q)^P+{ko,-k)YPL^+{qo,-q) • (14) 

Due to the form of the screening effects described in the previous section it is easy to see 
that the RG equations for the A-type couplings will not contain terms of order g"^. The Debye 
mass rriD in the longitudinal part of the gluon propagator (|^) removes any IR logarithm from 
the one gluon exchange diagram. 

In contrast, the B coupling may recieve such corrections. The leading contribution will 
come from the region of shell integration where one may neglect the Landau damping term 
G{Q) in gluon propagator (^. In this region the scattering angle satisfies the following 
condition: 1 — x^^^ < z = cos 6 < 1 — x'^/^, with x = i^—^- The angular momentum 
/-component Bi receives the contribution 

12 g^ 

3 3 4/12 7-1 l-z' 

(the factor of 1/3 corresponds to the attractive 3 channel, while the factor of 2/3 arises 
from the transversality of the magnetic gluon propagator). The correction is of the form 
6B1 ~ ^ proportional to the RG scaling parameter 6t = In Note that the leading 
term in 6B1 does not depend on the details of the angular cut-offs - that ambiguity is 
cancelled in the definition of t. 



The box diagram of figure ^ has the following integral representation^: 



A^) = / 7??i^A«P.-7''f(^/^ + P4)7'+p-7l7'^®7''k^/i-P4)7'-p-7l7% (16) 



^For notational simplicity, we suppress the chiral projectors Pl in our expressions, although they appear 
in our calculations. 
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where K^p^y^ incorporates both fermion and propagators (^: 



In order to evaluate the contribution to I(t) one should apply the familiar decomposition of 
momenta near FS: p = ftp + Pn, k = flk,q = flq with normal component being integrated in 
the limits A'^ < p,^ < . Note that due to screening effects the dominant regions of pn 
and p4 are around the origin, close to the FS. 

Consider, for example, the part of I{t) with both gluon lines being transverse. All other 
cases are done analogously. Capturing only the leading 0(ln^ A) behavior one gets 

ru\ 9^ f f dpidpn /• ,^ 7i7S-®7SV + li{p/ ■ ihj ® Y{p/ ' ih^ /1Q^ 
'^'^^ - WJ J^J'''^ ipl + pl)il - - ^2) ' ^''^ 

where zi is the cosine of the angle between flk and /i^, and Z2 is the cosine of the angle 
between /i^ and flp. Then, use the definition (^) of Bi in the expansion of transverse gluon 
propagators 

= -^E(2^ + l)^'^/W (19) 
1 - z 9 I 

on the interval^ 1 < 2; < 1 — x'^/^. Applying orthogonality conditions^ : 

j dn Pi (z,) Pv{z2) = 5u' Pi{z) 
I dn p,p, Pi {z,) Pv{z2) = 5iv ^-f Pi{z) , (20) 

one obtains the final answer for the /(t)T in the leading divergence approximation 

mr ^ -£^(5 74 ®7' + y7.®70^f(A)t ■ (21) 

The final answer for the (Minkowskian) RGE equations is 

dA^ ^ _Nf^^ 
dt 2 



- -(Af - 2AiBi + 5Bf 

dBl .n. o , .oo2\ , 9^ 



— (AJ - WAiBi + 13B?) + , (22) 

dt 6 V ' U 27/i2 ' ^ ^ 

where N = /i^/27r^. These RG equations are the same as obtained previously for the case 

of local interactions 0, except for the constant term in the Bi sector. A simple way to 

understand this is as follows: by expanding the four quark form factor in terms of orthogonal 

The longitudinal gluon propagators may be expanded in Legendre polynomials Pi{z) on the interval 
-1 < z < 1. 

"^The averaging over coordinate systems is assumed in the second orthogonality condition of (EQ). 
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local operators G;, we reduce the problem to one in which the gluon exchanges are replaced 
by local interactions. Hence the analysis of should apply and the same RG equations 
obtained for each value of 1. This is true as well for LR helicity scattering, which we do not 
consider explicitly here. Note that though the Ai couplings originate in the Debye screened 
sector, it would be a mistake to discard them in the RGE equations due to the fact that 
there is mixing with the Bi, which diverge near the FS. Diagonalization of ( p2D gives 



It is straightforward to solve these equations. For convenience, we define the spin combi- 
nation Si = Ai — 3Bi and the spin 1 combination Ti = Ai + Bi. We find that Si reaches a 
Landau pole at the scale 

Ar^f,g-'exp{-^) , (24) 

which is agrees with the result (0), despite differences between the RG equations of |^ and 
ours (see below). Ti, due to the opposite sign of the term, does not diverge and reaches 
the asymptotic value of Ti{t — > oo) = 

In order to compare our results with those of 0, it is necessary to convert four quark 
operators into scattering amplitudes. There are additional angular dependences introduced 
by the spin angular momentum of the quarks. Let us classify amplitudes by their total 
angular momentum, which is the sum of the spin and orbital components. Thus, for example, 
the j=0 channel receives contributions from both the 1=0, s=0 and 1=1, s=l operators. First 
we note that the spinor part of operators of type A introduce an additional factor of 2 cos^ | 
to the scattering amplitude, while the type B operators introduce a factor of 2 cos^ | + 
4 sin^ I . This leads to the following expression for the amplitude 

f{e) = J2{2i + i)[si + Ticose]Piicose) (25) 
I 



The components of total angular momentum j are given by 

— — 1 + — 



/. - + + I^T,,,] , (26) 



where we have used the identity 

{2l + l)zPi{z) = {l + l)Pi+i{z)+lPi^i{z) . (27) 

The RG equations for fj can be deduced easily from (^) and (^3]). For the lowest component, 
we have 

I/O ^ - ^ - A'So^ - . (28, 
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The equations for higher components are of the form 

dt^' ~ 27/^2 ^^^^ 3V2J + I ^-i + 2j + l ^-t-V • ^^^> 
Several remarks are in order. 

• These results are similar to, but different from, those of |]^. In particular, we do not find 
that the RG equations take the simple form [0 

= -^-Nff . (30) 

The rhs of our equation ( pOj) cannot be organized in terms of any simple amplitude fj. We 
believe that the treatment of quantum corrections in section 2 of that paper is too crude. 
In particular, iterating the amplitude f{p,k) in order to obtain the fj terms in the RG 
equation neglects some important spinor structure of the vertex which our calculation takes 
into account. 

• Fermi statistics constrain the operators as follows. We restrict ourselves to the color 
3 and isospin singlet channels, which are both antisymmetric. The remaining part of the 
wavefunction must be antisymmetric. In the antisymmetric s = channel, we must therefore 
have Z = 0, 2, 4, while in the symmetric s = 1 channel we must have / = 1, 3, 5, .... Thus 
the operators Si vanish for odd 1 and Ti vanish for even 1. From (^) we see that fj vanishes 
for all odd values of j. Note that this analysis is modified in the LR channels, which we do 
not consider here. 

• It is not necessary to solve the fj RG equations, as their behavior can be deduced from 
that of Si and T;. Since Ti never diverges, near the FS fj ~ Sj and diverges at the scale 

). Note that the ^Tf term cannot be neglected in (B^); its asymptotic value is — 



3 ^j ^c^x.x.vj^ ^.^f^.^^^^^ XX. y^^J, x^^ c^^jxxx^>.v^>.xv. y^.^^ x^ 27^3- 

• Because the Landau poles in all the j channels are the same, we naively expect to find con- 
densates with non-zero angular momentum, leading to the breaking of rotational symmetry 
[0. However, we will see in the next section by studying the gap equation that this is not 
the case. 



5 Rotational Symmetry Breaking and Gap Equation 

In this section we explore the issue of condensates with non-zero angular momentum using 
the gap equation. Our RG results suggest that the gap function could violate rotational 
invariance. Some recent papers have studied the s-wave condensate using the non-local 



gap equation, including the effects of magnetic gluons [16, 17, 18|. We will generalize their 
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approach to consider gap functions A(go, ^ which are functions of orientation. Remarkably, 
it is possible to show that in the leading order approximation only an s-wave condensate is 
formed. 

Let us introduce a two component field ^ = {ip, ip'^). The inverse quark propagator takes 
the form 

where A = 7oA"''7o. The gap is a matrix in color, isospin, and Dirac space, and in our 
analysis we allow it to depend on orientation. As discussed, the RG analysis shows that the 
condensate will form in the LL (RR) channels. Given this, the form of the gap matrix is 

ii 

= (A2)"'(r2).,C75 (A+(go, + a ■ q) + A_(go, (32) 

In our present weak coupling discussion, we are only interested in A+. A_ corresponds to 
a condensate of anti-quarks, and does not influence A_|_. Henceforth we shall only refer to 
A = A+. Making the usual FS approximations, the gap equation has the form 

. 2 f d^k D{qo- ko,q-k) A{ko,k,k) 
J (27r)4 e{k)^ + k^ + \A{ko,k,k)\^ 

where the interaction kernel D{qo — ko, q — k) is essentially the gluon propagator up to some 
additional factors arising from the particle projector and gamma matrices. Here k = \k\, 
g = |g| and e{k) = k — ft ; we assume ~ /^^ k q jji. For the next step in 

our analysis it is useful to separate the interaction kernel into angular momentum channels, 

D{q^~k^,q-k) = Y.{2l + l)Pi{cosek,) D\q^-k^,\q-k\) . (34) 

The coefficients are obtained by integrating the kernel against Legendre polynomials 
P;(cos^). Due to the divergence at small angle, all of the have the same value at leading 
logarithmic order. 

The gap equation can now be rewritten as 

A/ A -2^^1(0 A\ fJ^vi*fa ^\ -P'(go - ko, \q- k\) A{ko,k,k) 

im y^^r e^{k) + kl + \A{ko,k,k)\^ 

(35) 

For = -D(go — ^O; I"? ~ ^1) independent of 1, the sum over spherical harmonics reduces to a 
delta function: 6{(f)q — (pk) S {cos 9 g — cos9k). The angular integral on the rhs of the equation 
can be trivially evaluated, yielding 

. 2 f dkidkP D{qo- ko,\q- k\)A{ko,k,q) 
^^^°'^'^) = -^^ e^ik) + ki + \Aik,,k,q)\^ ■ ^''^ 
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This equation can be regarded as a set of identical equations, one for each orientation q. 
Since each equation is identical, the solution must be independent of orientation. 

This result could have been guessed directly from the fact that colinear scattering dom- 
inates the magnetic effects. In that approximation the kernel D is proportional to a delta 
function 5{k — g), and it is clear that (|33D only has an s-wave solution. 

Our analysis thus far has been within the leading logarithmic approximation. We can 
relax this condition by considering a gap function 

A(g)=A° + Ai(go,g) , (37) 

where A° is the constant (in orientation) solution obtained from (0), and A^ is a small 
perturbation which can depend on orientation. We will show that A^ is at most of order 
exp(— c/(7^) and hence negligible relative to the A°. Substituting (|37| ) into (|33|), we obtain two 
gap equations. The leading order equation is the just the usual one for an s-wave condensate, 
and determines A°. Note that we retain the complete interaction kernel -D(go — ko,q — k) 
here, without making the leading log approximation. The second equation contains terms 
of 0(Ai): 

Aif ^ 2 2 f dkodnkqD{qo-ko,ekq) \ |A°|2 1 

This equation was obtained after first performing the integral over k, in the approximation 
that A and D are slowly varying for k c:^ fi. The term in brackets in (|38|) suppresses the 
logarithmic divergence in the integral over k^ near the FS, although there is still a potential 
divergence from the small angle behavior of D{qQ — ko, Okq)- Hence the solution A^ is at most 
of order exp(— 1/(7^), and is negligible relative to A° in the weak coupling limit. 

The discussion in terms of A" and A^ is quite general. We can also apply it to lower 
densities, where the magnetic interaction is presumably screened by non-perturbative effects, 
and quark interactions are effectively local. In this case the largest condensate can be shown 
to be rotationally invariant |^, Because D in (^) is non-singular, a non-zero solution A^ 
only exists above some critical coupling (if at all). This is unlike the usual Cooper pairing 
instability in which an arbitrarily weak interaction can lead to a condensate. It seems that 
a large s-wave component tends to inhibit condensates of higher angular momentum. 

Having made some general observations about rotational invariance, we now concentrate 
on solving the gap equation explicitly, in order to check our result for the Landau pole (^). 
The authors of 0, |18[ do not consider higher orbital angular momentum components 
A;>o in their analyses, but it is straightforward to do so. 
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IS 



The explicit gap equation for Ai = A, neglecting the suppressed antiparticle contribution, 



A(a = f i-l^-g (oo) 

5 J {2nre-ik) + ki + \Aiko,k)\Hk-qr + G ' ^''^ 

Here we have written explicitly the factors resulting from the particle projector and gamma 
matrices. G is the Landau damping term which appears in the magnetic gluon propagator; 
we neglect the effect of the electric gluon. 

In order to consider higher angular modes of A we make the expansion: 

A{qo,q)=Y.{2l + l)Pi{q)Ai. , (40) 
I 

as well as a similar expansion of the gluon propagator with coefficients 
which are all of the same size at leading order. Using the identity 

o_ m=l 

PiicosO,,) = -— Y: YnekAkTYriO,,^,) , (42) 

ra=—l 

we obtain the following coupled gap equations 

2ig'^ f dkodk fi^ 



A, 



X 



Svr J (27r)2 e2(^) + ^2 + |^(^^^ ^)|2 
|aA, - ^^^^((/ + l)A+iA, + /A-iA,)] . (43) 



In obtaining this equation we have neglected the angular dependence of |A(fco,^)P which 
occurs in the denominator. This is justified if the gap turns out to be rotationally invariant, 
as expected from our previous arguments. It is easy to see that a self-consistent solution 
exists with all / > gaps zero, and the solution for A = Aq given by 



6 Discussion 



In this paper we investigated the renormalization group behavior of QCD at high density, 
concentrating on the effects of long range magnetic interactions. Our approach was some- 
what different from that of in that we focused on individual local operators rather than 
scattering amplitudes. The resulting RG equations are different, although the location of 
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the Landau pole is still given by (|l]). The disagreement results from two causes: (I) the 
transverse form of the propagator does not appear to have been used in leading to a 
different coefficient in the constant term of the RGE and (II) the treatment of spinor 
properties of the quantum corrections is different in the two calculations. We believe that 
the renormalization group applied directly to amplitudes does not properly compute the loop 
corrections. Some additional issues we attempted to clarify include the vahdity of the use of 
Landau damping in resummed gluon propagators, the gauge invariance of the computation 
and the size of subleading corrections. Our RG equations ( PBD are gauge invariant, and rep- 
resent the leading order result in a self-consistent expansion. Corrections to the coefficients 
in ( p3|) are of order 0{g) in the weak coupling limit. 

We also used the gap equation to investigate whether rotational symmetry is broken at 
asymptotically high densities. The gap equation analysis shows that scattering which is 
predominantly colinear leads to a rotationally invariant condensate. We found, in disagree- 
ment with a naive interpretation of the RG results, that any condensates of higher angular 
momentum are exponentially smaller than the s-wave condensate. We also checked that the 
solution of the gap equation agrees with the value of our Landau pole. 

The authors would like to thank Nick Evans, Deog-Ki Hong, James Hormuzdiar, Rob Pisarski 
and Dirk Rischke for useful discussions and comments. This work was supported in part 
under DOE contracts DE-FG02-91ER40676 and DE-FG06-85ER40224. 
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